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Abstract 

The model of an young star with a protoplanetary disk and a low-mass companion (q < 0.1), which 
is moving on a circular orbit, inclined to the disk plane, is considered. The hydrodynamic models of 
fs^ ' such a system were calculated by SPH method. It was shown the distortions in the disk, caused by the 



in 



orbital motion of the companion, lead to the strong dependence of illumination conditions of the disk on 

the azimuth (because of extinction between the star and the disk surface) and, therefore, it leads to the 

i-Q \ appearance of large-scale asymmetry in images of disks. The calculations showed the dependence of the 

»-^H. illumination on the azimuth was stronger in the central part of the disk than on the periphery. Bright and 

q ' dark domains are located not symmetry with respect to the line of nodes. The sizes and locations of the 

domains are depended on the model parameters as well as on the phase of the orbital period. The bright 

and dark domains do not follow the companion, but they make small amplitude oscillations with respect 

to some direction. The properties of the model, which were written above, open new opportunities of 

searching low-mass companions in the vicinity of young stars. The stars with protoplanetary disks, which 

are observed face-on or under the small inclination angle i, are the best ones for this purpose. 
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INTRODUCTION 

Many young stars are known to be surrounded by 
circumstellar gas-dust disks radiating in the infrared 
(IR) and submillimeter wavelength ranges (see, e.g., 
the review by Natta et al. (2000)). In the overwhelm- 
ing majority of cases, observations provide only in- 
formation about the integrated (i.e., from the entire 
disk) IR radiation (Mendoza 1966; Cohen 1975). In 
the visual spectral range, the circumstellar disks radi- 
ate very little, mainly through the scattered radiation 
of their central starqJ.Nevertheless, most of the cir- 
cumstellar disk images known to date have been ob- 
tained precisely in the optical spectral range with the 
Hubble Space Telescope by the method of coronagra- 
phy (McCaughrean and O'Dell 1996; McCaughrean 
et al. 2000; Grady et al. 2000; Krist et al. 2000; 
Clampin et al. 2003; and others). Among them, 
there are disks seen edge-on or at a small angle to 
the line of sight (Burrows et al. 1995; McCaughrean 
and O'Dell 1996; Stapelfeldt et al. 1998; and oth- 
ers) and disks seen nearly face-on (Grady et al. 2000; 
Krist et al. 2000; and others). Many of the disks 
have an asymmetric shape that is often explained by 
anisotropic scattering of stellar radiation by circum- 
stellar dust grains. This mechanism is efficient in the 
cases where the disk is inclined to the plane of the 
sky. 

However, the nature of the asymmetry in circum- 
stellar disk images is not exhausted by this mech- 
anism alone. The asymmetry of the disk surface 
brightness can also be caused by an anisotropic il- 
lumination of the disk by the central star, for exam- 
ple, due to the presence of spots on the star. This 
model was considered by Wood and Whitney (1998) 
in connection with the detection of variability in the 
image of HH30's protoplanetary disk by Burrows et 
al. (1996). It is quite real, because the rotational 
modulation of the brightness of young stars due to a 
nonuniform surface brightness is a well-known obser- 
vational fact (see, e.g., Vrba et al. 1986). 

Another cause of the anisotropic disk illumination 
can be the absorption of stellar radiation by circum- 
stellar matter if it is distributed nonuniformly in az- 
imuth. Hydrodynamic calculations show (Artymow- 
icz and Lubow 1996; Larwood and Papaloizou 1997; 
Sotnikova and Grining 2007; Hanawa et al. 2010; 
Kaigorodov et al. 2010) that such conditions can 
arise in young binary systems accreting matter from 
the remnants of a protostellar cloud (below referred 
to as a CB (circumbinary) disk). Under the action 
of periodic perturbations induced by the orbital mo- 
tion of the companion, spiral density waves and gas 
flows propagating into the system's inner regions are 



formed in the CB disk. If such a system is observed at 
a small angle to the disk plane, then, as our calcula- 
tions showed (Sotnikova and Grinin 2007; Demidova 
et al. 2010a, 2010b; Grinin et al. 2010), the motion 
of large-scale structures in it can cause great circum- 
stellar extinction variations and, as a consequence, a 
large-amplitude brightness modulation of the system. 
It turned out than a noticeable (in amplitude) photo- 
metric effect could be caused by the orbital motion of 
a low-mass companion or a protoplanet (Demidova et 
al. 2010b) and that this effect is enhanced if the com- 
panion's orbit is inclined to the disk plane (Grinin et 
al. 2010). 

It should be noted that the situation where the 
orbit of the companion (planet) does not coincide 
with the midplane of the common disk or the equato- 
rial plane of the central star is not something excep- 
tional. This can result from a nonaxisymmetric dis- 
tribution of specific angular momentum in the proto- 
stellar cloud from which the star and protoplanetary 
disk are born. A star whose rotation axis may not 
coincide with the rotation axis of the protoplanetary 
disk is formed during the gravitational contraction of 
such a cloud. This possibility has been discussed in 
recent years (see, e.g., Bate et al. 2010) in connection 
with the unexpected results of spectroscopic observa- 
tions of stars with exoplanets during the transit of a 
planet over the stellar disk. They showed that the 
orbital planes of several exoplanets do not coincide 
with the equatorial plane of the stars (Hebrard et al. 
2009; Winn et al. 2009a, 2009b; Pont et al. 2010; 
Narita et al. 2011; Brown et al. 2012). Previously, 
Burrows et al. (1995) found that the inner part of the 
circumstellar disk around /3 Pic Pic is inclined by sev- 
eral degrees to its outer part and hypothesized that 
this inclination was caused by the motion of a planet 
whose orbit is also inclined to the disk plane. The 
corresponding model was developed by Mouillet et 
al. (1997) and Larwood and Papaloizou (1997), and 
the planet itself has also been discovered recently (La- 
grange et al. 2009; Chauvin et al. 2012). The same 
picture is observed for the Herbig Ae star CQ Tau: 
the inner part of its circumstellar disk is inclined ap- 
proximately by 30° to the disk periphery (Eisner et 
al. 2004; Chapillon et al. 2008). 

In this paper, we investigate the influence of hy- 
drodynamic perturbations in such systems on the ex- 
tinction variations between the central star and disk 
surface. The goal of our calculations is to determine 
the disk illumination conditions at various distances 
from the center and its dependence on azimuth. 



1 UX Ori stars at eclipses (Grinin et al. 1991) and young 
objects with edge-on circumstellar disks constitute an excep- 
tion. 



THE METHOD OF CALCULA- 
TIONS 

We consider the model of a binary system that con- 
sists of a central star and a low-mass companion (g = 
to 2 : m\ < 0.1), moving in a circular orbit of radius 
a. The orbital plane is assumed to be inclined to the 
CB disk midplane. The mass of the star is taken to 
be 2 Mq this value roughly corresponds to the mass 
of a typical Herbig Ae star) ; the orbital period of the 
companion is 1 year. 

As in our previous papers, the 3D hydrodynamic 
calculations were performed by the SPH (Smoothed 
Particle Hydrodynamics) method in the isothermal 
approximation. The implementation of this method 
is described in detail in Sotnikova (1996). One of the 
main hydrodynamic parameters, the effective viscos- 
ity of the disk, was determined via the dimensionless 
speed of sound c, expressed in fractions of the Keple- 
rian velocity in the companion's orbit. This param- 
eter was varied within the range c = 0.001 — 0.08. 
The parameters of the problem also include the mass 
accretion rate onto the binary components M a and 
the opacity k, per gram of circumstellar matter (as 
in our previous papers, k was taken to be 250 cm 2 /g 
(Natta and Whitney 2000), which corresponds to the 
optical properties of the circumstellar dust in young 
circumstellar disks in the visual spectral range) . The 
dust-to-gas mass ratio was taken to be the same as 
that, on average, in the interstellar medium, 1 : 100. 

The calculations were performed using from 60000 
to 2 x 10 5 test particles in a region of radius 10 a. Our 
calculations showed that after 30 rotations the sys- 
tem finishes the relaxation stage (with the formation 
of a central gap) and the distribution of SPH parti- 
cles reflects best the behavior of the matter in the 
CB disk. The disk model calculated in this way was 
smoothed over the cells of a 3D mesh (with a step of 
0.1a). This procedure allowed the influence of ran- 
dom fluctuations in the distribution of SPH particles 
to be reduced while retaining all details of the flow 
structure. 

As an example, Fig. 1 shows the distribution of 
SPH particles in one of the calculated models (be- 
fore the smoothing procedure described above). We 
see that the deviations from azimuthal symmetry are 
greatest near the inner boundary of the CB disk (the 
disk is puffed up). At distances of the order of five or- 
bital radii, the distribution of matter in the CB disk 
is already azimuthally symmetric, and such a disk dif- 
fers little in structure from the disk around a single 
star of the same mass. To investigate the influence 
of density waves and nonuniform structure of the in- 
ner region of the CB disk on the illumination of its 
outer layers, the model distribution of particles was 
smoothly extended to distances of 50 a by using the 



standard model of a flared disk: 



h(r)=h ( 



r y 
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where 2 h(r) — is the geometrical thickness of the disk 
at distance r from the symmetry axis, the parameter 
j3 = 5/4 (Kenyon and Hartmann 1987), and the ini- 
tial values of ho = a, and the radius ro = 5a. 



Calculation of the Disk Illumination 

We calculated the disk surface illumination from the 
formula 

F(r ^ ) = 4^FT7?) sin7 ' (2) 

where L* — is the luminosity of the central star, 
r(r, <fi) — is the optical depth of the dust layer be- 
tween the star and an arbitrary point on the disk at 
distance r from the symmetry axis and azimuth 4>, 
7 — is the angle between the vector connecting the 
star and the point on the disk surface and the tan- 
gent to the disk surface at this point (for details, see 
Tambovtseva et al. 2006). 

When calculating r(r, 4>) we used the same method 
as in our previous papers (see, e.g.,Demidova et al. 
2010a): having compared the accretion rate of test 
particles onto the binary components with the ac- 
cretion rate specified as a parameter of the problem 
(in the range M a = 10~ 9 - lO~ 11 M yr~ 1 ), we de- 
termine the test particle mass m<j. Then, r(r, <fi) — 
rrid ■ n{r, (p) ■ k/S, where n(r, </>) is the number density 
of test particles in the column between the star and 
a point on the disk surface with coordinates r h <j>, S 
is the cross section of this column. 



Allowance for the Scattered Radiation 

The direct stellar radiation is the main disk surface il- 
lumination source. This radiation does not penetrate 
into the shadow zone because of strong absorption on 
the way, in the perturbed disk region, and it is illu- 
minated only by the radiation scattered by circum- 
stellar dust. To estimate the contribution from this 
radiation, we calculated the intensity of the scattered 
radiation for several models in the single-scattering 
approximation. This approximation is applied in the 
cases (see, e.g., Tambovtseva et al. 2006) where the 
single-scattering albedo differs noticeably from unity, 
because under these conditions the role of multiple 
scatterings is minor. 

Below, we used the same scheme of calculations 
as in the above paper (Tambovtseva et al. 2006): 
the inner CB disk region (r < 5a) was divided into, 
cells, each being considered as a scattering element 
with density pi = m<j ■ n%, where n, s the number of 



particles in the zth cell. For each surface element of 
the CB disk, the scattered light was calculated as the 
sum of the fluxes scattered by each cell toward the 
chosen disk surface element. We took into account 
the absorption of radiation by dust both on the way 
from the star to the scattering volume and on the 
way from it to the illuminated disk surface element: 



F sc (r,4>) = k s 



Li* 



/ — J Andj And 



/(^)e- T *cos&. (3) 



Here, d\ is the distance from the primary compo- 
nent to the scattering cell, di is the distance from 
the scattering cell to the unit CB disk surface ele- 
ment illuminated by the scattered radiation, k sc is 
the scattering coefficient, t^ = Ti + t^ is the sum of 
the optical depths of dust on the way from the pri- 
mary component to the scattering cell (n) and from 
the scattering cell to the surface element on the outer 
CB disk (Y2), £j is the angle between the direction of 
the scattered radiation from the center of the scatter- 
ing cell to the surface element on the CB disk under 
consideration and the normal to the disk surface at 
this point. In our calculations, we used the Henyey- 
Greenstein phase function: 



/(*) 



1 



1-.9 2 



(4) 



4tt(1+ 5 2 - 2. 9 cos 0) 3 / 2 ' 

where 9 is the scattering angle, g is the asymmetry 
factor (taken to be 0.5). The single-scattering albedo 
was taken to be 0.4, which roughly corresponds to the 
optical properties of circumstellar dust in the V band 
(Natta and Whitney 2000). Taking this into account, 
the scattering coefficient is k sc = 100cm 2 /g. 



RESULTS 

We calculated the disk surface illumination by the 
method described above for several models. The in- 
put model parameters are the orbital inclination to 
the disk plane a, the component mass ratio q, and 
the disk viscosity parameter c. Figure 2 presents the 
results of our calculations for two models differing by 
the orbital inclination: a = 10° and a = 30°. We see 
that, in both cases, the disk illumination is highly 
nonuniform in azimuth due to the absorption of the 
direct stellar radiation in the perturbed region (near 
the inner disk boundary) . The degree of deformation 
of the inner CB disk increases with a causing the 
area of the shadow from it on the disk to increase. 
A nontrivial peculiarity in the disk brightness distri- 
bution attributable to the distribution of matter in 
the perturbed CB disk region is that the line sepa- 
rating the bright and dark regions on the disk does 
not coincide with the line of nodes. For quantitative 



estimations, we determined the location of this line 
in such a way that the ratio of the integrated illumi- 
nations of the disk regions on different sides of this 
line was maximal. Our calculations showed that the 
inclination of this line to the line of nodes changes 
with model parameters within a comparatively nar- 
row range: <j>i w 50 — 70°. In particular, in the models 
with a = 10° the angle 4>i = 57°, with a = 30° - 
02 = 59° (Fig. 2). 

As can be seen from Fig. 3, the location of the 
shadow zone relative to the line of nodes in the mod- 
els considered is almost independent of the compan- 
ion's orbital inclination. The illumination minimum 
is reached in a region with an azimuthal angle 4> s h w 
330°. Since the minimum illumination is due to the 
maximum absorption of radiation in the perturbed 
disk region, the result obtained suggests that the re- 
gion of maximum disk deformation does not coin- 
cide with the region in which the companion rises 
to the maximum height above the CB disk plane 
(4> = 270°). Such a distribution of matter in the inner 
disk stems from the fact that the secondary compo- 
nent entrains part of the matter once it has passed 
the point of greatest distance from the CB disk plane, 
which causes the CB disk thickness in this direction 
to increase. Other things being equal, the larger the 
viscosity the higher the efficiency of this process (see 
below). 

The deformation of the inner CB disk region also 
depends on the companion's mass. The greater the 
value of q, the larger the amplitude of the perturba- 
tions in the disk and the stronger the asymmetry in 
its illumination. In this case, the position of the il- 
lumination minimum still corresponds to cj> s h ps 330° 
(Fig. 4). 

Viscosity affects the illumination conditions in the 
following way: the CB disk becomes denser and thin- 
ner with decreasing viscosity. As a result, the area of 
the shadow on the disk is reduced. For example, in 
the models with a "cold" disk (c < 0.03) the disk is il- 
luminated uniformly already at R > 20a, while in the 
models with c c = 0.05 (at the same value of q) the 
shadow zone extends to R = 30a. The location of the 
line separating the bright and dark regions on the disk 
also depends on viscosity and varies within the range 
4>i = 50 — 65°. In addition, the azimuthal angle of the 
illumination minimum in the shadow zone increases 
with viscosity (Fig. 5). This suggests that the disk 
matter is entrained more actively by the secondary 
component during its orbital motion and shields the 
star more strongly with increasing viscosity. 

Figures 6 and 7 show the two limiting situations 
considered in our calculations where a noticeable asym- 
metry in the CB disk illumination is still possible: 
the model with a small orbital inclination to the disk 
plane (a = 3°) and the model in which the mass of 
the perturbing body is of the order of several Jupiter 



masses (q — 0.003). Under these conditions, the per- 
turbations in the disk induced by the companion's or- 
bital motion are small. Nevertheless, as can be seen 
from Figs. 6 and 7, an asymmetry in the disk illumi- 
nation due to an asymmetric distribution of absorb- 
ing matter is clearly seen in the central disk region 
in these cases as well. Our calculations showed that 
the degree of asymmetry in the disk illumination is 
more sensitive to variations in a, than to variations 
in the companion's mass. Therefore, even a planet 
with a mass of several Jupiter masses can produce a 
noticeable asymmetry in the protoplanetary disk il- 
lumination if its orbit has an appreciable inclination 
to the disk plane (a > 10°). 



rotation^ (Wood and Whitney f998). In our case, as 
a result of the companion's orbital motion, the bright 
and dark regions execute small synchronous oscilla- 
tions on the disk surface, with their shape chang- 
ing only slightly. This is clearly seen from Fig. 10, 
which shows two disk models. The companion is at 
an orbital phase of 90° , in one model and at the op- 
posite end of the orbit, at a phase of 270°, in the 
other model. We see that the boundary between the 
bright and dark zones in these two limiting situa- 
tions changed only slightly. The position angle of the 
line separating these two regions changes with orbital 
phase, though within a narrow range (Fig. 11). 



Influence of the Scattered Light 

In addition to the direct radiation, the disk surface is 
also illuminated by the scattered light. The contribu- 
tion from this radiation on the bright (illuminated) 
part of the disk is negligible compared to the direct 
stellar radiation. However, the scattered light in the 
shadow region makes a tangible contribution to the 
disk surface brightness. Where, according to Figs. 3- 
5, the illumination by the direction radiation is close 
to zero, the disk surface brightness is nonzero due to 
the scattered light. 

We calculated the disk illumination in the single 
scattering approximation for several models. Two of 
them are presented in Fig. 8. Analysis of the calcu- 
lated models showed that the scattered light reduces 
the depth of the shadow zone only slightly in the case 
of small disk deformation (at small a). However, at 
a = 30° the contribution from the scattered light in 
the shadow zone is more significant. As can be seen 
from Figs. 8 and 9, a bright spot emerges in the 
shadow zone in the inner CB disk. It appears, be- 
cause the inner CB disk is strongly deformed and has 
a significant thickness at large angles a. Therefore, 
the scattered light from the perturbed region on the 
CB disk illuminates the outer part of its surface at a 
less acute angle than in the models with a lower a. As 
a result, in the models with strong disk deformation, 
the scattered light makes a tangible contribution to 
the illumination of the shadow zone. 



Orbital Motion of the Companion and 
the Disk Illumination 

The changes in disk illumination conditions with or- 
bital phase are of special interest. Calculations showed 
that these changes occur by no means in the same way 
as, for example, in the model of a spotted star, where 
the regions of maximum and minimum illumination 
on the disk move azimuthally, following the stellar 



Asymmetry of Disks During Corona- 
graphic Observations 

As has been pointed out in the Introduction, the 
method of coronagraphy, whereby the central star 
(along with the central region on the disk) is cov- 
ered, is commonly used in the observations of circum- 
stellar disks. We modeled this situation. Figure 12 
shows two disk illumination distributions in the same 
model. In one of them, the central region of the sys- 
tem is covered by an opaque shield. The shield size 
relative to the disk size was chosen in such a way that 
it roughly corresponded to the conditions of Hubble 
Space Telescope observations of the circumstellar disk 
around one of the nearest Herbig Ae stars, HD 163296 
(see Wisniewski et al. 2008). We see that, in this 
case, the shadow zone is covered almost completely by 
the artificial shield. Consequently, the coronagraphic 
method of imaging widely used in studying circum- 
stellar disks can hide valuable information about the 
existence of an asymmetry on the disk caused by a 
low-mass companion from the observer. 

It should be noted that in close binary systems like 
those we consider here with an orbital radius of the 
order of several AU, the dust in the perturbed region 
of the CB disk near its inner boundary makes a major 
contribution to the circumstellar extinction. Other 
potentially important extinction sources include the 
circumstellar disk surrounding the central stars and 
the gas flows in the inner region of the system. In the 
implementation of the SPH method (with a constant 
smoothing length) we used, this region is modeled 
rather roughly. Since this region has small sizes in 
close binary systems like those considered here, its 
contribution to the extinction is negligible. In wider 
pairs, the inner region of the system can contain much 
dust, and more accurate methods should be used for 
its calculations (Hanawa et al. 2010; Fateeva et al. 



2 The shadow region created by a dusty disk wind from 
the secondary component in the models by Tambovtseva et 
al. (2006) moves over the disk in the same way. 



2011; Sytov et al. 2011; de Val-Borro et al. 2011). 



CONCLUSIONS 

Our calculations show that the perturbations in the 
circumstellar disk of a young star induced by the or- 
bital motion of a low-mass companion whose orbit 
is inclined to the disk plane can cause a noticeable 
asymmetry in the disk illumination. It differs fun- 
damentally in properties from what is given by the 
previously considered models. For example, in the 
models of disks with anisotropic scattering of stellar 
radiation by dust grains, a disk asymmetry is possible 
only when the disk is inclined to the plane of the sky 
(Augereau et al. 1999). During face-on observations, 
such a disk will be azimuthally uniform. In contrast, 
in our case, an asymmetry in the disk illumination 
exists at any disk inclination and is best seen during 
face-on observations. 

In contrast to the models in which the moving 
bright and dark regions on the disk follow the rota- 
tion of a spotted star (Wood and Whitney 1998) or 
the motion of a companion that is the source of a 
disk wind (Tambovtseva et al. 2006) , the orbital mo- 
tion of the companion in our model does not cause 
any synchronous motion of the shadow on the disk 
but leads only to small oscillations of the bright and 
dark regions.. To a first approximation, the boundary 
between them retains its location in projection onto 
the plane of the sky. In this case, the line separating 
the illuminated region and the shadow zone, coincide 
with the line of nodes. The position angle of this line 
measured from the line of nodes is 40-60° (depending 
on the model parameters). 

The amplitude of the perturbations in the disk in- 
duced by the orbital motion of the companion, along 
with the parameters of the bright and dark regions 
on the disk, depend on the component mass ratio, 
viscosity, and orbital inclination. Other things be- 
ing equal, the smaller these parameters, the smaller 
the shadow zone on the disk. The minimum mass of 
the companion at which the asymmetry in the disk 
illumination is still noticeable in our calculations is 
about six Jupiter masses (Fig. 6). 

Allowance for the influence of scattered light showed 
that it reduces the depth of the shadow on the disk 
without changing the asymmetry in the disk illumina- 
tion. The case where the companion's orbit is highly 
inclined to the disk plane (Fig. 8) constitutes an ex- 
ception. In this case, the scattered radiation forms a 
bright spot in the shadow region due to strong defor- 
mation of the inner CB disk. This effect was obtained 
in the single-scattering approximation and deserves 
a more detailed study by solving the 3D radiative 
transfer problem. 



The degree of asymmetry in the illuminated re- 
gion on the disk is maximal in the central part of the 
disk with a radius of 20-30 orbital radii and decreases 
with increasing distance from the center. During 
coronagraphic observations, this region is covered by 
a mask (an artificial moon) covering the star. We ex- 
pect a big breakthrough in studying the asymmetry 
of circumstellar disk images once the ALMA interfer- 
ometer has been put into operation. It will operate 
in the submillimeter wavelength range and will allow 
the central regions of protoplanetary disks to be seen 
without any distortions. 
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Figure 1 : Distribution of matter in the disk after 30 rotations of the system: (a) face-on view, (b) and (c) 
disk sections along the y and x axes; here, a is the angle between the orbital plane of the binary system 
and the CB disk plane at the initial instant of time. The distances are in units of the binary's semimajor 
axis. The model parameters are: the component mass ratio q = 0.1, the viscosity parameter c = 0.05 and 
a = 10°. The angle <j> is measured from the y axis (coincident with the line of nodes) in the direction of the 
companion's rotation. 
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Figure 2: Top: illumination of the outer part of the CB disk by the binary's primary component for two 
models differing by the orbital inclination of the companion to the CB disk plane. The model parameters 
are: q = 0.1, c = 0.05, a = 10° (left) and a = 30°. (right). The readings along the x and y axes are in units 
of the binary's semimajor axis; the color scale of illumination is in erg/(sxcm 2 ). Bottom: the same in the 
form of isolines; the straight line corresponds to the line separating the bright and dark (shadow) regions. 
The dotted lines in the lower graphs indicate the regions where the illumination is greater than that of the 
outer disk at a distance of 50 a; the solid lines indicate the regions where it is smaller than this value. 
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Figure 3: Normalized dependences of the disk illumination (F n ) on azimuth at three distances from the 
center (for the same models as in Fig. 2). The distance R is in units of the semimajor axis, a — 30° (top) 
and a = 10° (bottom). 
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Figure 4: Same as Fig. 3 for three models differing by the companion mass: q = 0.01 (top), q = 0.03 
(middle), and q = 0.1 (bottom). In all three models, c = 0.05, a = 5°. 
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Figure 5: Same as Fig. 3 for three models differing by the viscosity: the upper, middle, and lower plots 
correspond to c = 0.05, c = 0.04 and c = 0.02, respectively. In all three models, q = 0.1, a = 5°. 
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Figure 6: Small perturbations in the limiting cases. Top: the disk illumination for the models with q = 0.003, 
c = 0.05, a = 10° (left) and q = 0.1, c = 0.05, a = 3° (right). Bottom: the same in the form of isolines; the 
inclination of the line separating the dark and bright regions is <j>i = 46° (left) and 4>i — 44° (right). 
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Figure 7: Same as Fig. 3 for two models with small perturbations: the upper and lower plots correspond 
to the models with q = 0.003, c = 0.05, a = 10° and q = 0.1, c = 0.05, a = 3°, respectively 
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Figure 8: Top: disk illumination including the scattered light for the models with q — 0.1, c — 0.05, a = 5° 
(left) and a — 30°. (right). Bottom: the same in the form of isolines; the inclination of the line separating 
the dark and bright regions is cf>i = 49° (left) and (f>i — 59° (right). 
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Figure 9: Same as Fig. 3 for the models with and without the scattered radiation. The model parameters 
are: q — 0.1, c = 0.05, a = 5° (the upper pair of plots) and a = 30°. (the lower pair of plots). The upper 
and lower plots in each pair show the dependence F n without and with the scattered light, respectively. 
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Figure 10: Illumination of the outer CB disk (with the scattered light) for two positions of the secondary 
components in its orbit. Left: the companion is at = 90° (under the disk), the inclination of the separating 
line is (j>i = 58°; right: <j> = 270° (the companion is above the disk), 4>i — 71°. The model parameters are 
q = 0.1, c = 0.05 enda = 10°. 
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Figure 11: Inclination (</>;) of the line separating the bright and dark regions versus orbital phase. The 
model parameters are q = 0.1, c = 0.05 and a = 10°. Phase corresponds to the position of the component 

at <p = 0°. 
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Figure 12: Illumination of the outer CB disk. The model parameters are q — 0.1, c = 0.05 and a = 10°. 
The entire disk and only its outer region (the inner part with a size R < 20a is covered by an opaque shield) 
are shown on the left and the right, respectively. 
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